WHAT IS CLAIMED IS: 

1 . A statistical facial feature extraction method, comprising: 

a first procedure for creating a statistical face shape model based on 
a plurality of training face images, including: 
5 an image selecting step, to select N training face images; 

a feature labeling step, to respectively label feature points located in 
n different blocks of the training face images to define corresponding shape 
vectors of the training face images; 

an aligning step, to align each shape vector with a reference shape 
10 vector to thus obtain aligned shape vectors; and 

a statistical face shape model computing step, to use a principal 
component analysis (PCA) process to compute a plurality of principal 
components based on the aligned shape vectors to create a statistical face 
shape model, wherein the statistical face shape model represents the shape 
15 vectors by combining a plurality of projection coefficients; 
and 

a second procedure for extracting a plurality of facial features from 
a test face image, including: 

a test face image selecting step, to select a test face image; 
20 an initial guessing step, to guess initial positions of n test feature 

points located in the test face image, wherein the initial position of each test 
feature point is a mean value of the feature points of the aligned shape 
vectors; 

a search range defining step, to define n search ranges in the test 



face image, based on the initial position of each test feature point, wherein 
each search range corresponds to a different block; 

a candidate feature point labeling step, to label a plurality of 
candidate feature points for each search range; 
5 a test shape vector forming step, to do combination of the candidate 

feature points in different search ranges in order to form a plurality of test 
shape vectors; and 

a determining step, to match the test shape vectors respectively to 
both the mean value and the principal components for computing a 
10 similarity, and to accordingly assign one feature point corresponding to one, 
having the best similarity, of the test shape vectors as facial features of the 
test face image. 

2. The method as claimed in claim 1, wherein in the feature labeling 
step of the first procedure, the feature points are coordinates for comers of 

15 eyes and mouth on each training face image. 

3. The method as claimed in claim 1, wherein the feature labeling 
step of the first procedure manually labels the feature points of each training 
face image. 

4. The method as claimed in claim 1, wherein the reference shape 
20 vector is one of the shape vectors. 

5. The method as claimed in claim 1, wherein the aligning step of 
the first procedure uses a 2D scaled rigid transform algorithm to align each 
shape vector with the reference shape vector. 

6. The method as claimed in claim 5, wherein the aligning step of 



the first procedure further comprises the steps of: 

selecting the reference shape vector as x i = (x iI ,y il ,...,x in ,y iB ) and 
one of the shape vectors as x } = (x^ , y jf x jn , y jn ) ; 

computing a squared Euclidean distance E between the vectors Xj 
and xj based on the following equation 
E = (Xj - M (N) (a,0)[ Xj ] - t) T ( Xi - M^Ca^XxJ - t) , where 

M (N) (a,0)[x.]- t is a geometric transform function defining with a 
plurality of transfer parameters to align the vector x J5 the transfer parameters 
include a rotating angle 9 , a scaling factor a , and a shifting vector 
10 represented by t = (t x ,t y ), and M (N) (a,0) is a 2nx2n diagonal blocked 

^occos0x jk -asin0y jk ^ 
v asin0x jk +acos0y jk y 



matrix as well as M(a,0) 



M(a,0) = 



L y jk . 
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for 1 ^k^n 5 as 



finding the smallest squared Euclidean distance and corresponding 
rotating angle Q } , scaling factor and shifting vector represented by 
15 tj =(t xj ,t yj ) to align the shape vector xj so similar as the reference shape 
vector Xj; 

computing a sum of smallest squared Euclidean distances after the 
N shape vectors are all aligned so similar as the reference shape vector, 
ending the aligning step when the sum is smaller than a predetermined 
20 threshold; 

computing a mean value of the feature points in each block for the 
aligned shape vectors to define a mean shape vector for each aligned shape 
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1 N 

vector as x = — Z x a > wherein x a is the aligned shape vector; and 

assigning the mean shape vector as the reference shape vector and 
the aligned shape vectors as the shape vectors and then repeating the 
aligning step until all shape vectors are aligned. 
5 7. The method as claimed in claim 6, wherein the transfer 

parameters is obtained by a least square algorithm. 

8. The method as claimed in claim 1, wherein the statistical face 
shape model is a point distribution model (PDM). 

9. The method as claimed in claim 8 5 wherein the statistical face 
10 shape model computing step of the first procedure further comprises the 

steps of: 

computing a mean value of the feature points of the aligned shape 

-In 

vectors to define a mean shape vector as x = — £x a , wherein x a is the 

N a=l 

aligned shape vector; 
1 5 subtracting each aligned shape vector by the mean shape vector to 

form a matrix as A = |d v ,d_ ,...,d xvr I, wherein d x = x - x ; 

computing a covariance matrix of the matrix A as C = AA T ; and 
computing a plurality of principal components according to 
eigenvectors derived from Cv s k = X\v s k with eigenvalues corresponding to 
20 the covariance matrix C, to form the statistical face shape model, wherein 
X\ represents eigenvalues of the covariance matrix C, v s k represents 
eigenvectors of the covariance matrix C, and 1 ^k^m, where m is the 
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dimension of the covariance matrix C for X s > X s > ... > X s . 

12 m 

10. The method as claimed in claim 1, wherein each shape vector xj 
consists of n feature vectors Sj located in different blocks, so an average 

1 n 

value as t = — of the feature vectors Sj corresponding to special blocks 

5 of all shape vectors Xj is defined as a feature template. 

11. The method as claimed in claim 1, wherein in the initial 
guessing step of the second procedure, scaling of initial guess shapes 
formed by the test feature points is aligned so similar as the test face image. 

12. The method as claimed in claim 10, wherein the candidate 
10 feature point labeling step of the second procedure further comprises the 

steps of: 

ic 

labeling a plurality of reference points derived from I; = t + £bp. 

j=i J J 

respectively in each search range, where t is the feature template of block 
corresponding to a search range, pj is j-th principal component of the 
1 5 statistical face shape model computed from the training feature vectors, and 
bj is an associated projection coefficient; 

k 

using e =|| Ij ~ t - ZbjP; ll 2 to compute an error value between one 

of the reference points and the corresponding principal component pj and 
projection coefficient bj; and 
20 selecting preceding k smallest error values and defining the k 

smallest error values as feature points of the search range. 



13. The method as claimed in claim 12, wherein the test shape 
vector forming step of the second procedure does combination of the 
candidate feature points in different search ranges to thus form k n test 
shape vectors. 

5 14. The method as claimed in claim 10, wherein the determining 

step of the second procedure further comprises the steps of: 

using the average value of the aligned shape vectors and the 
principal components of the statistical face shape model to represent an 

k 

approximate value of the test shape vector as x = x + Zb?p x , where x is a 

1 0 mean shape vector defined according to the mean value of the feature points 
of the aligned shape vectors, p* is j-th principal component of the statistical 
face shape model, and b* is a corresponding projection coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape 



vector represented by x = M(a,G) 



x + Zblp 



+ t , where 0 , a and t are a 



1 5 rotating angle, a scaling factor and a shifting vector respectively; 

computing a normalized distance of the aligned test shape vectors 



by d(x) = Jz 



'bl v 



; and 



assigning one candidate feature point corresponding to one, having 
the smallest normalized distance, of the aligned test shape vectors as facial 
20 feature of the test face image. 

15. The method as claimed in claim 10, wherein the determining 
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step of the second procedure further comprises the steps of: 

using the average value of the aligned shape vectors and the 
principal components of the statistical face shape model to represent an 



— k 

approximate value of the test shape vector as x = x -f £bj p x , where x is a 

j=i 1 



5 mean shape vector defined according to the mean value of the feature points 
of the aligned shape vectors, p* is j-th principal component of the statistical 
face shape model, and bj is a corresponding projection coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape 



vector represented by x = M(a 3 0) x + Zbjp- J + t , where 0 , a and t 



are a 



10 rotating angle, a scaling factor and a shifting vector respectively; 

computing an error value between the test shape vector and the 
average value of the aligned test shape vectors by 

e(x) = w, i|| I, (x) - t, - tb)p] || 2 +w 2 d(x) , where £|| I. (x) - 1, - ±b)p] ll 2 . 

is a similarity of the test shape vector to corresponding aligned shape vector 
15 x a , and d(x) is the normalized distance of the aligned test shape vectors; and 
assigning one candidate feature point corresponding to one, having 
the smallest error value, of the test shape vectors as facial feature of the test 
face image. 

16. The method as claimed in claim 15, wherein the error value is 



20 computed by an equation e(x) = w, 



v v 
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+ w 2 d(x), 



